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ABSTRACT 


The following is a report on the design and construction of a 
fixed-base variable-stability simulator facility combined with the task 
of landing to a carrier or runway. 

The solution was mechanized on a hybrid computer with the analog 
computer solving the equations of motion and the digital computer used 
for storage, control and graphics generation. The display was in the 
form of a computer-drawn picture on a graphics terminal. Control was by 
a simulated cockpit placed in front of the display and connected to the 
analog computer. 

Dynamic validation was considered excellent with the modal periods 
of the simulated aircraft agreeing very closely with those of the actual 
aircraft. 

The visual display was deemed yery good as sufficient visual cues 
were provided to enable consistent landings by experienced pilots. 

This project was undertaken not as a design of training aid but 
rather as a research tool for further studies in control systems, human 


engineering and aircraft dynamics. 
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І. INTRODUCTION 


The purpose of this project was to construct a multiple task simu- 
lator, an inherent quality of which had to be the ease of its extension 
aS a research tool to various real world problems that arise. In addi- 
tion, the use of interactive computer graphics would be studied as a 
means of visual display for the simulator. The decision was made to 
design the simulator to the initial task of approach and landing to an 
aircraft carrier deck. | 

The hub of the simulator system is a hybrid computer, part digital, 
part analog. The analog computer solved the equations of motion while 
the digital computer was used for control, storage and graphic build-up. 
The display was handled by a graphics processor and terminal connected 
directly to the digital computer. Control of the simulator was through 
a Simulated cockpit placed in front of the graphic display, with stick 
and throttle outputs tied directly to the analog. This then formed the 
control loop as diagrammed in Figure l. 

At the beginning, a decision on the allocation of equipment func- 
tions had to be made. Several all-digital moving-base simulators were 
already constructed [Ref. 1]. As space, money and time were limited, 
moving base was out of the question but a pure digital solution was not. 
On the other hand the department already had a "Mini Link" Simulator 
utilizing an all analog solution. However, due to the lack of non-linear 
equipment and the need for generating a visual perspective display, the 
pure analog solution could not be implemented. Then too, the digital 
solution was too slow for the needed integration and did not provide 


for a pilot-station computer interface. The obvious solution was the 
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Figure 1. Lines of Communications. 





hybrid computer. The non-linearities of the aircraft coefficients, the 
resolution of the forces and moments and the generation of the display 
would be tasked to the digital computer. The integration of the equa- 
tions of motion and the interfacing would be done by the analog. In this 
manner the best features of both computers could be utilized. 

The display would be a computer-drawn picture on a graphics terminal 
connected to the digital computer. Construction of the picture would be 
done by the digital computer wnile the display and refresh would be done 
by the graphic digital processor. 

The "cockpit" was simply a seat and attached controls placed in 
front of the graphics CRT. Stick and throttle outputs from the cockpit 
were then fed directly into the analog computer, thus closing the 
control loop. 

The design features incorporated into the Carrier Approach Landing 
Simulator (CALS) were: 

1. Ability to change the aircraft being simulated 

2. Capability of simulating non-linear aerodynamic data 
Automatic scaling of the analog computer solution 
Fixed-base inside-out display 


Ability to expand or change the simulator task 


a ο A O 


Simplicity of setup and operation 
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IT. THE SIMULATOR 


The simulator is capable of accurately simulating vehicle dynamic 
response in conjunction with a dynamic inside-out display with its view- 
ing plane affixed to the vehicle. The analog solution is updated at a 
rate of 18 samples per second which leads to excellent dynamic response. 
A new graphic display image is generated at the rate of 18 frames per 
second with a digital processor refresh rate of 40 frames per second. 
This leads to a display with virtually no flicker and in general very 
good graphic qualities. 

The simulator was designed to accept the following inputs in the 
form of punched cards: the coordinates of the runway; the runway cons- 
tants (glide slope and lens placement); the aircraft constants (eight, 
moment of inertia, wing area etc.); initial conditions of the problem 
(initial position, altitude and velocity) the maximum expected range 
of the state and control variables in the form of scale factors: and the 
aerodynamic coefficients of the simulated aircraft. The preparation of 
this data deck is outlined in Appendix F. 

The simulator then processes this data, autoscales the analog 
computer based on the given and self-generated scale factors, sets the 
potentiometers of the analog and places the initial conditions on the 
proper integrators. Control of the simulator is then transferred to 
the pilot who responds as per the displayed instructions. At the 
completion of a run the simulator analyzes the landing and displays the 
results. This analysis includes sink rate at touch down, line-up, 
airspeed and type of landing (crash, bolter or arrested landing on a 
certain wire). At this time options are offered to either fly again 


with the same parameters, change the parameters or stop the program. 
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Besides these options, the operator can suppress printed output if 
desired, select a yaw damper (ON/OFF) and inhibit translational movement 
while flying the aircraft. Also, at the disposal of the operator, is ап 
eight-channel strip recorder and an oscilloscope that are patched into 
the analog computer. This gives the operator the ability to output the 
time history of any variable he wishes. 

The cockpit consists of a chair upon which the throttle and control 
stick have been mounted. The stick, mounted on the right arm, is a 
Gemini control stick providing all three controls: yaw, pitch and roll. 
At the top of the stick is located one of the control buttons used in 
the control of the program from the pilot station. Mounted on the left 
arm is a throttle plate containing a small throttle and another control 
button. The outputs of the cockpit are fed via a wire bundle to the 
rear of the analog computer for processing. 


A complete operating manual is contained in Appendix H. 
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ПІ. SIX DEGREE OF FREEDOM AIRFRAME EQUATIONS 


R.M. Howe [Refs. 2 and 3] discusses three coordinate systems - body, 
wind, and stability axes - and concludes that the use of wind axes for the 
translational equations and of body axes for the rotational equations of 
motion is the best choice. As can be seen from Figure 2, the body axes 
differ from the stability axes by the aircraft angle of attack, while 
the stability axes differ from the wind or flight path axes by the side 
slip angle. 

Assuming that the three angular rates of the aircraft (P,Q,R) 
about the body axis are given, the anaular rates about the stability 


axis (Р.,0.,К.) сап be computed (Eqn 101-103). 


p. Picos ca + IR Sin; o (101) 
Q. = Q (102) 
R. = P sina + Q cos a (103) 


These angular rates can then be normalized using equations 104-106. 


` 


P_b Ä 
(2 ТЕЛЕ (104) 
S 2V 
Q b 
E. 5 
Qs 7 2V (105) 
R b 
"αι - ---- (106) 
ς 2ү 


The aerodynamic forces expressed along the stability axes can then 


be calculated (Eqn 107-109). 


F 
X (x ) 
aero _ D 
^m Чүл (107) 








AL - Body Axis 


CAPS - Stability Axis 


AL - Wind or Flight Path Axis 


Figure 2. Aircraft Coordinate Axis. 
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The forces in the body axes are then resolved into forces in the 
stability axes (Eqn 110-112). Thrust and the force of gravity are 


included at this time. Thrust was assumed to have only an X component. 


Fy Fy 
st ; : aero 
ET С UE g sin ο) ον απ ος 6 005 πο en (110) 
Fy F 
S : aero 
ы ος ο δη τσ: (111) 
Ε F 
7 Z 
— = -(Т. g sin a) sin αν ο σος ο ος о كق ۽‎ (112) 


The forces are then resolved from the stability axes to the wind or 


flight path axes (Eqn 113-115). 





i n i | 

ΠΠ π᾽ ο05 πο ο (113) 

Mw X sin 8 * Ys cos 8 (114) 
m m m 

F. F 
Z Z 

ES - 5 (115) 
m m 


From these equations the translational equations of motion can be 


determined. Solving for the derivatives yields equations 116-118. 


F 
X 
V Е W (116) 
m 
ET : 
Z P sin β 
ES. м4. 5 м ΓΕ, 
mV cos β cos 8 S 








T. W 
Sr р (118) 


Further, the Euler angular rates can then be calculated (Eqn 119- 


121). 
. А 
-р = - (R er > sin $) (119) 
-0 =- Q cos ç + R sin à (120) 
-p Py sin 0 (121) 


Now focusing our attention on the moment equations, the three equa- 


tions about the stability axes can be written as follows (Eqn 122-124): 


P_SbC R _SbC daSbC 6rSbC 
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M. 5-6, ποῦς. m е 
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These moments are then resolved into the body axes yielding 


Equations 125-127. 


L L. N. И 

ет 05 а - р sino (125) 
XX XX 27 

M 

τ---τ-- (126) 
yy yy 

N Ls і ПЕ 

E = r sin α + 7 Cos а (127) 
22 XX 72 


The rotational equations of motion can then be written (Equations 


128-130). 
Г, „(А + РО) 


! (I, - 1.) 
dE уу a E Al 
ХХ XX XX 


- 


(128) 
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As suggested by Howe [Ref. 2], a simplified method for computing 


the velocities in the inertial frame yields equations 131-133. This 


approach uses small angle approximations in its development. 


The 
can then 
exist at 

P = 

(| = 


С 
— 


V cos 8 Cos y (131) 
V cos € sin y | j 
V(-sin 0 + sin a cos 8 cos ¢) (133) 


various derivatives that were solved for in the above equations 
be integrated yielding the state variables that were assumed to 


the start (Eqns 134-145) 


SP dt (134) 
sQ dt (135) 
JR dt (136) 
fa dt (137) 
f8 dt (138) 
fy dt (139) 
je dt (140) 
sg dt (141) 
SV dt (142) 
/S, dt (143) 
S, dt (144) 
rS, dt (145) 





The equations themselves, when examined, seemingly divide themselves 
into two groups. One that can be solved by digital means and one that 
can be solved by the use of the analog computer. In fact this logical 


division is the one chosen for the simulator. 
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IV. GRAPHICS PRESENTATION 


Since the primary scan during a carrier approach consists of only 
three things, the "meatball" or the Fresnel Lens Optical Landing System, 
the aircraft's angle of attack and the line-up with the runway center- 
line, it was decided to incorporate only these three visual cues into 
the display. 

The angle of attack was represented by an angle of attack indexer 
situated in the upper right hand corner of the display. Its operation 
is exactly the same as the one found in the jets/turboprops operating in 
the Navy today. The indexer has five descrete states as shown in 


Figure 3. 





States ] 


Figure 3. Angle of Attack Indexer States. 
Optimum/Landing angle of attack was represented by State 3. As one 
increases the angle of attack of the aircraft, State 2 would appear 
indicating the aircraft was a "little slow" and eventually State 1 would 
appear indicating that the aircraft was "slow". Proceeding the other 
direction from State 3, that is decreasing the angle of attack, first 
State 4 would appear indicating the aircraft was a "little fast" and 


finally State 5 would appear indicating the aircraft was "fast". In 
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terms of degrees, assuming optimum angle of attack to be 0°, the 


following would exist: 


STATES RANGE OF AOA 
1 Greater than 2.0? 
2 Between 1.0° and 2.0° 
3 EO που ος 
1 Between -1.0° and -2.0° 
У Less than -2.0° 


The Fresnel Lens is represented by two fixed rectangular boxes for 


the datum lights and a movable square for the "meatball" as shown in 


REE: 


Figure 4. Centered Meatball. 


Figure 4. 





The operation is again the same as found in the Navy systems. If 


the "meatball is above the datum bars, as in Figure 5, the aircraft is 


porem 


Figure 5. High Meatball. 


above the glide slope. 





If the "meatball" is below the datum bars, as in Figure 6, the 


aircraft is below the glide slope. 


— 


Figure 6. Low Meatball. 
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The object is to keep the "meatball" centered within the datums. 
The "meatball" has a 1.5 degree visibility arc centered on the optimum 
glide slope and will disappear if this is exceeded. Also, once you fly 
past the lens itself the "meatball" will disappear. However, in this 
simulation the datums will remain visible at all times, only the 
"meatball" disappears. 

Since the lens is affixed to the carrier/runway in real life, the 
size of the presentation will change. However, since the vast majority 
of flight cues are absent in the simulator, the size of the simulated 
Tens will remain fixed. This presentation is placed in the upper left 
hand corner of the display. 

The line-up cue is received from the picture of the runway itself. 
The picture shown is as if a movie camera were affixed to the nose 
(X-axis) of the aircraft. The field of view is a square with an 
angular limit of plus or minus 18.5 degrees in the vertical and 
horizontal directions. Within this window lies the horizon and the 
runway outline. Tic marks at 0° and plus and minus 10° of pitch are 
given along the left side of the window. As the aircraft maneuvers 
through space the runway and horizon move dynamically as if filmed 


by the movie camera. The complete display is shown in Figure 7. 
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Figure 7. 





The Complete Graphics Display. 
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Ve PROCESSING THE DISPEAY 


The processing of the angle of attack indexer is straightforward as 
it is solely a function of the angle of attack of the aircraft. A simple 
test of the angle of attack will yield one of the five previously 
mentioned states. 

The Fresnel Lens or "meatball" is again a simple task to process. 
Knowing the x-coordinate and the z-coordinate of the aircraft and the 


x-coordinate of the lens, the glide slope angle can be calculated. 


(5.26) 






i) 
Figure 8. Glide Slope Computation. 
The angle GS in Figure 8 is the glide slope angle while Up £p! and 


(X, 2, ) are the coordinates of the plane and lens respectively. 


65 - tan (2,-2,)/06-X)) 


Angle GS is then calculated and compared to the optimum glide slope, 
and the "meatball" is drawn accordingly. If GS is greater than the 
allowable limits, or if Хр « Xy > the "meatball" is not displayed. 

The processing of the runway is by far the most complicated. How 


to transform a three-dimensional object like the runway into a two- 


dimensional object capable of being displayed on the screen of a CRI? 
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The answer was found in two works, one by L.G. Roberts [Ref. 3] and the 
other by R.B. Desens [Ref. 4]. Roberts introduces the idea of a homo- 
geneous coordinate system which introduces a fourth coordinate or 
variable scale factor. He also constructs a single transformation 
matrix for all points. This transformation matrix (Ех) те а 
product of five matrices: rotation, translation, perspective, another 
translation and scale. Each point is post-multiplied by the H-Matrix 
and then the scale factor is divided out yielding the display coordinates. 
Desens used this to process pictures of an aircraft carrier from various 
positions. He also reduced the H-Matrix to a product of three individual 


matrices as shown in Figure 9. 


А, А ДА АА, АЛ 1 ЕСЕ 1 -S/F 





za D 3 
y B. В, В; 0 BB, BB, 88. 0 1 0 0 
С, % C4 0 CC4 CC» CC. 0 0 ] 0 
Nez dl 0 0 0 (aea Te 5 | 
Object Axis Viewing Plane Perspective 
Rotation Orientation Scale 
(Euler Angles) Offset 
Viewing Plane 
Translation 


Figure 9. H-Matrix 
Figure 10 has the two coordinate systems, that of the earth and 
that of the aircraft. The viewing plane or camera is placed on the 
nose of the aircraft while the runway is oriented in the object or 
earth reference system. The first matrix is made up in part by the 
direction cosines of the Euler angle rotation of the object (runway). 
In this case the runway is fixed. The Euler angle direction cosines 


are given below and are based on the rotations shown in Figure 11. 
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Figure 10. 


Graphic Coordinate System. 


Euler Angle Direction Cosines 


- COS ß cos 
= cos g sin 


= -sin ß 


а 


а 


- -COS y Sin a +sin y sin В соѕ а 


= COS a COS y +sin y sin 8 Sin o 
- Sin y cos 8 
= sin y sin a teos y Sin B COS o 
--Sin y COS a *COS y Sin 8 Sin a 
= cos B COS y 
2 
α 
Order of Rotation Y 
αν Вот Y 
B 
X 
Figure 11. Euler Angle Rotation. 
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Xr Үт апа Z+ are the translation distances between the two axis 


systems. 


m^ object ¬ view plane 
n object ¬ vie plane 
object y view plane 

The second matrix is made up entirely of the viewing-plane-orienta- 
tion direction cosines. The direction cosines are given below in terms 
of the Euler angles of the aircraft. 

AA, = COS y Cos 0 

AA, = sin y coS 9 -cos y sin 8 sin à 

AA, = sin y sin à *cos y sin 8 cos q 
BB, = -COS 8 sin y 
BB, = cos y cos ф +sin 8 sin ψ sin ¢ 


BB, = cos y Sin ọ -sin y sin 0 cos q 


CC, = -sin 0 
CC, = -COS 8 sin & Order of Rotation 
CC, = cos 0 cos à y (yaw), ə (pitch), e (roll) 


The third matrix consists of the offset option, La and Yo? which 
allows one to look at an offset portion of the viewing plane. This 
option is not used on this display. The matrix also contains the scale 
factor (S) and the focal length (F). The scale factor is set at 1/2 
and the focal length, the distance between the viewer and the viewing 
plane, is set at 1.5 ft. The field of view is controlled by S and F 
and is equal to сап (S/F). 

A point on the object (X,Y,Z) is converted to homogeneous coordinates 
by the addition of another scale factor (W) yielding (X', Y', Z',W) 


where X' = WX, Y' = WY, Z' = WZ. This then is post multiplied by the 
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H-Matrix yielding (X'', Y'', Z'', W). The display coordinates of the 


processed point can then be found by dividing by W' 


Mtr! = ІІ! 
т! = αυ 
The coordinate X''' - X''/W' is an indicator of depth and can be 


used for depth cues, however, it is not used in this display. Each 
point of the runway is passed through the H-Matrix yielding the display 
coordinates of each point in turn. 

Turning to the problem of the horizon and referencing Figure 12 


below the following computations are presented. 


A 





Figure 12. Horizon and Earth. 


For an aircraft at point A with altitude z above the earth with 


radius r_, the true horizon is at point B. This is depressed a distance 


e? 
Zz, from the plane DE which is tangent to the earth at point D. However, 
0 
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a pseudo horizon exists at point E which always has a z coordinate of 
zero and visually coincides with the true horizon but has the advantage 
of being quicker to draw. 


From Figure 12, since triangle ADE and ABC are similar this yields: 





- oz) 
AB 
(z + © = AB“ tr 2 


1 
(2 - i 


assuming er, to be very much greater than z yields: 


— 


(41.8)2 3 


1 
a = = K(z)2 where K == 0 
ee) 





Nj = 


The depression angle of the horizon is: 
dep = 90° - tan la) = 90° - tan" (K(2)2/2) 

Referring to Figure 13, which is a view of the window and horizon, 
the object is to calculate YMID, the y intercept of the horizon. The 
distance PITCH en be calculated knowing the pitch angle and the field 
view. 

PITCH 7» o/Field of View 

TPITCH is PITCH minus the depression angle DEPA. Simple trig yields 

cos 4 = TPITCH/YMID, where 4 is the roll angle of the aircraft. 


YMID = TPITCH/cos o 
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Knowing YMID and the roll angle the horizon can then be constructed. 
Upon the completion of the runway and the horizon, the lines are placed 
through a software window algorithm that will cause only the lines with- 


in the window/field of view to be displayed. This computation is done 


for every iteration of the problem solution. 
| 
| 
| 


1 з SS 


da α 
кї 
ғ 


YMID 





Figure 13. Calculating Horizon Intercept. 
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VI. THE PROGRAM 


As stated in Section III, a logical division of tasks exist. To 
integrate the derivatives of the twelve state variables by digital 
algorithms takes too much computational time, therefore the logical 
thing to do was to D/A these variables and integrate them on a real- 
time basis using the analog computer. The results are then A/D'd back 
to the digital computer for calculation of the derivatives. This forms 
the computational loop. | 

When a real-time simulation is undertaken the time involved for the 
solution of the problem must be a minimum. Great pains were taken to 
insure that this was indeed done. The following is a tabulation of the 
significant time-saving steps that were involved: 

1. Resolution of the angles in the digital computer using a 

sine table lookup. 

2. Fast linear interpolation used in both the sine table 

lookup and the aerodynamic coefficient lookup. 

3. Preprocessing of the aerodynamic data to yield usable 

aerodynamic coefficients. 

4. General and liberal use of program constants to combine 

two or more arithmetic operations. 

5. Use of the SQRT function vice the exponential function. 

6. Multiple assignment. 

7.  Rewriting by the Laboratory staff of the A/D - D/a routine 

resulting in the SUBROUTINE ADDA. 

8.  Rewriting of the digital processor control program 


GATED to yield a faster version GATDI. 
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9. The use of equivalencing scalars to vector elements to 

reduce the indexing time within the digital computer. 
10. Judicious programming and a dedicated cleaning up process 
of the generated code for the computational loop. 

The total time for the computational loop usina the above procedures 
was cut from an original 190 milliseconds to 55 milliseconds. This gives 
an-approximate solution update rate and graphic frame rate of 18 samples 
or frames per second. 

The simulator was finally flown but the dynamic response was not 
quite as expected. The aircraft seemed to be unstable in the Dutch Roll 
mode, so a yaw damper was constructed. This helped, but at times the 
aircraft was still unstable. Investigation into dynamic error analysis 
[Ref. 5] demonstrated that the computational loop time delay was causing 
a degradation in the damping ratio. The lowest damping ratio was in the 
Dutch Roll mode, hence the appearance of the instability there. The 
solution was to predict ahead a value to be used by the program using 
the known derivatives. When this prediction and update was incorporated, 
all instabilities were corrected. 

Engine dynamic response was added in an effort to make the throttle 
response realistic [Ref. 6]. Also, in an attempt to increase the number 
of cues available to the pilot, variable frequency noise was generated 
in proportion to the power applied. When weighed against the results 
achieved, it was decided to discard the noise generator in favor of the 
Simplicity of setup and operation. 

The digital program in its entirety can be found in Appendix A, 
while the analog diagram can be found in Appendix B. In an effort to 


make the program more readable, Appendix C contains descriptors of all 
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the program і 
program variables used, and Appendix D contains a program flow 


diagram. 
B Sample output of the program can be found in Appendix G 
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VIT. CONCLUSIONS 


Dynamic validation of the simulator was done by comparing the various 


modes of response. The tabulation below shows that in fact the dynamic 


response of the simulator is excellent when compared with the actual 


Bera ft. in this case the A-7. 


σι ὃν C гә 


Simulator Aircraft 
Dutch Roll Natural Frequency .25 cps .25 cps 
Short Period 2.0 sec 2.3 sec 
Phygoid Period 37 sec 38 sec 
Roll Mode Heavily damped in both cases 
Spiral Divergence Time Constant 29 sec 25-35 sec 


The simulator has now met the self-imposed requirements set forth 


in the Introduction. Briefly reviewing the points in the Introduction: 


i. 


2. 


E 


Multi-task simulator - Currently three follow-on studies are 
being conducted using the simulator. These include a spin 
simulation, investigation into stall dynamics and a study 
involved with a LSO work station. 

Ability to change the aircraft being simulated - Demonstrations 
have been conducted with the "standard" Navy runway and the 
carrier runway dimensions as well as linearized and non- 
linearized aerodynamic data. Time involved to change, in each 
Case, was approximately one minute. 

The use of interactive graphics for the display - Very good 
results were achieved using the inside-out display. Results 


show that there is a learning curve as with any simulator but 
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sufficient visual cues are provided to enable consistent landings 
by experienced pilots. 

4. Simplicity of setup and operation - The number of dials or 
buttons have been kept to a minimum; as a result the system can 
be up and operating within five minutes. As for the operation, 
the simulator has been flown by all manner of people with little 
or no explanation needed. 

Possible areas of follow-on study might include the following: 

1. The study of Heads Up Displays (HUD's), since only the graphic 

presentation need be changed. 

2. Autopilot/Gust Alleviators/Stability Augmentor systems that 

can be modeled on an analog board can be incorporated, as approximately 
one-half the board 15 left unused for this purpose. 

3. Force-stick analyzation with a realistic tracking problem. 

4. Re-design of the Optical Landing System to reduce the false 


motion cues that currently exist in the present system. 
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APPENDIX A - THE DIGITAL PROGRAM 


This appendix includes listings of the main program, written in 
FORTRAN, and the SUBROUTINE ADDA, the A/D - D/A subprogram. The sub- 
routine is written in METASYMBOL, the assembly language for the 


XDS - 9300 and was prepared gy the computer laboratory staff. 
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APPENDIX B - THE ANALOG PROGRAM 
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APPENDIX C - PROGRAM VARIABLES 


The following is an alphabetic listing of the program variables used 

in the simulator program. 

A - Array of angles, in order, alpha, beta, psi, theta, and phi 

A10 - Initial conditions on A(1), the angle of attack 

AÀ - Array of analog to digital variables 

AD - Array of analog to digital variables 

ADOTB - Normalized time derivitive of the angle of attack 

ADOTN - Array of the negatives of the time derivitatives of the angles 
in A 

ADOTO - Temporary storage for ADOTN during the predict and update phase 
of the program 

ALI - Aerodynamic rolling moment about the body axis 

ALPHA - Real argument based on angle of attack and used for coefficient 
lookup 

ALSI - Aerodynamic rolling moment about the stability axis 

AMI - Aerodynamic pitching moment about the body axis 

ANI - Aerodynamic yawing moment about the body axis 

ANSI - Aerodynamic yawing moment about the stability axis 

AOA - Actual angle of attack in degrees 

AOAF - Angle of attack break point for fast cheveron 

AOALDG - Optimum landing angle of attack 

AOALF - Angle of attack break point for fast cheveron and donut 

AOALS - Angle of attack break point for slow cheveron 

ARG - Argument used during sine and cosine lookup (type real) 


ASZ - Absolute value of SZ 
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B - Wing span 

B2 - One half the wing span B 

B3 - Scale factor divisor used to divide out the scale factor of the 

homogeneous transformation coordinates 

BETA - Real argument based on side slip angle used in coefficient 
lookup 

C - Array of aerodynamic coefficients for a given angle of attack and 

sideslip angle 

CA - Array of cosines of the angle array A 

CB - Mean aerodynamic chord 

CB2 - One half the mean aerodynamic chord CB 

CCA(n) - Scalar equivalents of the array CA 

CHOYC - Choice of program options offered by main program 

COEFA - Array of coefficients dependent on angle of attack only 

COEFAB - Array of coefficients dependent on angle of attack and side 

slip angle 

CONI - Program constant equal to the number of degrees in a radian 

CON2 - Program constant equal to .2 divided by .75 

CON3 - Program constant used in computing the depression angle of the 
horizon 

CON4 - Program constant equal to the sine of 30 degrees 

CON5 - Program constant equal to the cosine of 30 degrees 

CON6 - Program constant equal to the ratio of CON3 to CON5 

D - Constant equal to 2.5 times the sine of the roll angle and used in 

calculating the horizon 
DA - Aileron deflection 


DEGA - Degrees of angle of attack for which COEFA and COEFAB are tabulated 
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DEGB - Degrees of side slip for which COEFAB are tabulated 


DEGMAX - Max angle of the visibility cone formed by the square viewing 


window 


DIT - Flying tail deflection 


DITO - Initial flying tail deflection setting 


DR - Rudder deflection 


F - Focal length of viewing image 


FXAM - 


FXSM - 


FYAM - 


FYSM - 


FYWM - 


FZAM - 


FZWM - 


Aerodynamic force in the X-direction divided by the mass of the 
aircraft 

Force in X-direction of the Stability axis divided by the mass 
of the aircraft 

Aerodynamic force in the Y-direction divided by the mass of the 
aircraft 

Force in the Y-direction of the stability axis divided by the 
mass of the aircraft 

Force in the Y-direction of the wind axis divided by the mass 
of the aircraft 

Aerodynamic force in the Z-direction divided by the mass of the 
aircraft 

Force in the Z-direction of the wind axis divided by the mass 


of the aircraft 


G - Acceleration due to gravity 


GEE - Scaled acceleration due to gravity 


GS - Glide slope in degrees 


GSO - Optimal glide slope 


H - Graphic transformation matrix (H-Matrix) or array 


HI - First component array of the H-Matrix 
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H2 - Second component array of the H-Matrix 

H3 - Third component array of the H-Matrix 

HHnn - Scalar equivalent of the H-Matrix or array 

HINSCT - Horizontal intersection of a line, used in the software window 

HMIR - Height of Lens below the landing surface 

I - Integer variable usually a counter 

IA - Do-loop counter based on angle of attack 

IAGN - A counter to check the number of times a point has gone through 
the software window 

IALPHA - Integer conversion of ALPHA 

IARG - Integer conversion of ARG 

IB - Do-loop counter based on the side slip angle 

IBETA - Integer conversion of BETA 

ICHOYC - Integer conversion of CHOYC 

IDE - Graphic array containing the dynamic portion of the display 

IDEV - Graphic device number (1 or 2) 

IDIR - Graphic block directory for the graphics digital processor 

IDON - Array of packed points for the donut portion of the indexer 

IER - Error flag returned by graphic subroutines 

IFAST - Packed graphic array of the fast cheveron 

IPLS - Integer counter equal to I plus one 

ISF7 - Integer scale factor based on thrust divided by mass 

ISLO - Packed graphic array of the slow cheveron 

ISQ - Array/graphics block of fixed data, basically the square or window 

ITDIR - Text directory for the graphics digital processor 

ITEXnn - Line of text that is sent from the digital computer to the 


graphic digital processor 
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INIRE - Wire caught on arrested landing 

IX - Operating values of start/end points (X-coordinate) used in 
software window 

IX1 - Equivalent of XSTART, used in software window 

IX2 - Equivalent of XEND, used in software window 

IY - Operating values of start/end points (Y-coordinate), used in 

` software 

IY1 - Equivalent of YSTART, used in software window 

IY2 - Equivalent of YEND , used in software window 

J - Do-loop counter | 

KK - Array of integers used in outputting coefficients 

KPOT - Array of integers used in outputing analog pot settings 

LL - Landing result flag 

NAD - Number of analog to digital conversion variables/trunk lines 

ND - Display flag based on present condition of the angle of attack 
indexer 

NDA - Number of digital to analog conversion variables/trunk lines 

NULL - Octal value used to null out a line of text 

P - Angular velocity of the aircraft about X-axis 

Pn - Array of runway points 

PDOTN - The negative of the time derivitive of P 

POT - Array of pot settings for the analog 

PPnn - Scalar equivalents of Pn 

PS - P about the stability axis 

PSB - PS normalized 

Q - Angular velocity of the aircraft about the Y-axis 


QDOTN - The negative of the time derivitive of Q 
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05 - 0 about the stability axis 

QSB - QS normalized 

QUE - Dynamic pressure 

R - Angular velocity about the Z-axis 

RAMP - X-coordinate of the ramp of the runway 

RDOTN - The negative of the time derivitative of R 

RHO - Density of air 

RIXX - Moment of inertia of the aircraft about the X-axis 

RIXZ - Product of inertia about the X and Z-axis 

RIYY - Moment of inertia of the aircraft about the Z-axis 

RMASS - Mass of the aircraft 

RS - R about the stability axis 

RSB - RS normalized 

S - Wing area of the aircraft 

SA - Array of sines of the angle array А 

SDOT - Array of time derivitives (velocities) of the distances in all 
three inertial directions 

SFl - Scale factor based on maximum velocity in X-direction 

SF2 - Scale factor based on maximum velocity in Y-direction 

SF3 - Scale factor based on maximum velocity in Z-direction 

SF4 - Scale factor based on maximum angular velocity in X-direction 

SF5 - Scale factor based on maximum angular velocity in Y-direction 

SF6 - Scale factor based on the maximum angular velocity in the Z-direction 

SF7 - Scale factor based on thrust divided by mass of the aircraft 

SF8 - Scale factor based on maximum angle of attack 

SF9 - Scale factor based on the maximum side slip angle 


SFIO - Scale factor based on the maximum yaw angle 
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БЕЗ8 - Scale factor based on ratio of SF5 to SFII 

SF39 - Scale factor based on ratio of SF4 to SFl2 

ο οὐ - Ratio of SFZ to SFI 

SFA - Array of scale factors, made up to SF8 through SF12 
SFACTR - Scale factor of the H-Matrix 

SLOPE - Slope of a display line, used in software window 

SNKRT - Sink rate of the aircraft in ft/sec at the conclusion of a run 
SSAn - Scalar equivalents of SA 

SX - Actual distance of aircraft in X-direction (inertial) 

SXO - Initial condition on SX 

SXN - Nromalized Distance of aircraft in X-direction (inertial) 
SY - Actual distance of aircraft in Y-direction (inertial) 

SYN - Normalized distance of aircraft in Y-direction (inertial) 
SZ - Actual distance of aircraft in Z-direction (inertial) 

SZO - Initial condition on SZ 

SZN - Final cutoff value on SZ (height of eye) 

T - Thrust of the aircraft 

TALPHA - Tangent of the yaw angle 

TAU - Time delay of the prediction and update 

TM К Thrust divided by mass of aircraft 

TPn - Temporary point number n 

TRIG - Array of sines of various angles 

TX - X-direction value used to chop the runway 

V - Velocity 

VO - Initial condition on velocity V 

VDOT - Time derivitive of V 


VDOTO - VDOT of the iteration before 


91| 





VINSCT - Vertical intersection of a line, used in window 

W - Weight of the aircraft 

Wn - Wire position on runway 

WIND - Velocity of the wind in knots 

X - X value used in software window 

XEND - Array of X-coordinate ending points of display lines 

XMIR - X position of the lens 

XSTART - Array of X-coordinate starting points of display lines 

XTEMP - Temporary storage for X 

Y - Y value used in software window 

YEND - Array of Y-coordinate ending points of display lines 

YMB - Y-coordinate of the "meatball" 

YMID - Y intercept of the horizon 

YO - H-Matrix offset in Y-direction 

YSTART - Array of Y-coordinate starting points of display lines 

YTEMP - Temporary storage for Y-coordinate in software window 

Z - Temporary storage array used in multiplying H1, H2, and H2 to yield 
the H-Matrix 


Z0 - H-Matrix offset in Z-direction 
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APPENDIX D - DIGITAL PROGRAM FLOW DIAGRAM 











RUNWAY 
DATA 


ESTABLISH 
PROGRAM 
CONSTANTS 


BUILD FIXED 
PORTION OF 
GRAPHICS 





J3 






SCALE 
FACTORS 






| RUNWAY 
| CONSTANTS 







COMPUTE 
NEW SCALE 
FACTORS 





SCALE 
FACTORS 









AIRCRAFT 
COEFFICIENT 


М, 2 ι 


AIRCRAFT 
COT FICIENT 











COMPUTE USABLE 
AIRCRAFT 
COEFFICIENTS 







USABLE 
AIRCRAFT | 
COEFFICIENTS § 


COMPUTE 
POT 
SETTINGS 








94 








SINE AND 
COSINE 
LOOKUP 






ANALOG ` 
IN COMPUTE? 





COBREICIENT 
LOOKUP 
ANALOG TO 
COMPUTE 









RESOLVE 
AND FORCES 
MOMENTS 









INITIALIZE 
AND CHOP 
RUNWAY 






| PREDICT 








UPDATE 


95 





BUILD 
TRANSFORMATION 
MATRIX 


BUILD ARRAY 
μπε» 


BUILD 
HORIZON 







SOFTWARE 
WINDOW 





BUILD 
MEATBALL 








96 


BUILD 


RUNWAY AND 
HORIZON 
DISPLAY 













BUILD 
ANGLE ОБ 
ATTACK 
DISPLAY 


COMPLETE 
"GRAPHICS 
_ DISPLAY 






BUILD 
LANDING 
RESULTS 











LANDING 
RESULTS 


FLY AGAIN 








BUILD 
PROGRAM 
OPTIONS 







CHANGE 
PROGRAM 
VARIABLES 


97 


APPENDIX E - OPERATING MANUAL 


The program exists on tape in two forms, one is in a Source deck 
format (SI) and must be compiled prior to execution. This takes approxi- 
mately 15 minutes. This tape is stored in the tape locker and is titled 
CALS NEW VERSION. The other version is a core dump (SAVE) version. Just 
prior to the execution of the program, after the compilation and loading 
phases are complete, a copy of core is dumped on tape. Execution of 
this version (RERUN) is approximately 10 seconds. Again this is stored 
in the tape locker under the title of CALS CORE DUMP. 

All the required cards for the program are filed in the file under 
CARRIER APPROACH LANDING SIMULATOR. The cards filed include the master 
version of the program, several data decks and the control cards for 
running any of the different programs. The control cards are color 
coded with explanations of their use written on each set. 

The instructions for operating the digital computer and graphics 
display are already published and can be obtained from the Laboratory 
Office. 

The following is a sequence of events and instructions for operating 
the simulator assuming all equipment is turned on. 

1. Load the version of tape selected on the tape drives. Select 

the proper mount (2 for SI, 3 for RERUN). 

2. Load and ready the card reader with control cards followed 

by the DATA. 
3. Ready the Line Printer. 
4. On the Digital Main Frame, Punch IDLE, RESET, RUN, CARDS. 


Cards and tape should start being read at this time. 
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Load the Analog Program board and the Analog Logic board 
into the computer (Boards number 7). 
On the Analog Keyboard punch KEYBOARD, LOCAL, POTSET. In the 
far right control box, only the following can be lit: IDCXI or 
IDCX.1 and REAL TIME. All other lights must be punched out. 
When complete punch DIGITAL CMPTR. 
On the Analog the following options exist: 
DSO - UP - no printed output 

CENTER - printed output 
р51 - ОР - distance integrators disabled 

CENTER - distance integrators enabled 
DS2 - UP - yaw damper ON 

CENTER - yaw damper CFF 
DOWN positions on these switches are momentary contact positions. 
On the Analog: 
Turn the inner dials of DFl2 and DF13 clockwise until resistance 
ìs felt. Then set the outer dial of ОРЕ12 їо .1 апа the cuter 
dial of DF13 to 10. 
Under Master Clock (Analog) punch RUN. The run light on the 
light panel should be on and the DF13 light should be flashing. 
At the Graphics terminal teletype to be used (usually number 2), 
load the controling program by typing 

RESET ("GATDI", 120) ' 
If the above is answered with 
FILE NOT FOUND 

type БЕЗБЕН | “САТЫ 20)! 
Then execute the program by typing GATED: At this point 


nothing should be on the display screen. 
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12. 
19. 
14. 


Що, 


18. 


19. 


20. 


21. 


By this time the Digital computer should have compiled and 
loaded the main program and the following will be typed out on 
the Digital consol teletype: 

INPUT AGT NUMBER 
The input light will then come on. To respond, type the AGT 
number (1 or 2) of the unit being used, followed by a carriage 
return. 
The Title should appear on the selected AGT. 
The Data will be read in from the card reader. 
If output is selected, the output will be printed on the line 
printer. 
The Analog Pots will be set by the Digital computer. The 
Address and Ratiometer readings should be changing. 
Set the cockpit in front of the AGT (will only reach AGT-2). 
The cockpit is stored behind the Analog computer. 
The Instructions should flash on the screen. 
Punch the button on the Throttle plate and the display will 
appear. 
Punch the button on the control stick to fly. If during a run 
you wish to abort that run, punch the button on the Throttle 
plate. 
At the completion of the run the Landing Results will be 
displayed, followed by a short delay in which to read them. 
At this time two options are offered. To fly again punch the 
button on the Control stick and the program jumps back to the 


Instructions (#17). 
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С: 


28. 


To receive the expanded program options, punch the Throttle 
button. Follow the displayed instructions and type your 
selection on the AGT teletype unit. This will result in one 
of three things happening. 

1. Logical STOP to the program. 

2. Jump to the Instructions (#17). 

3. Jump to the selection of the AGR number (#11). 
At the completion of the computer time. Return all used 
equipment (cards, tape, boards, cockpit) to their storage 


locations. 
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APPENDIX F - PREPARATION OF DATA DECK 


The following is a guide to be used in preparing the data deck to 
be used with the simulator. 

Cards 1 - 6 - The first six data cards consist of the actual 
coordinates of the six runway points. The first ten columns of the card 
are the x coordinate in feet, the second ten columns are the y coordin- 
ate in feet, the third ten columns are the z coordinate in feet and the 
next ten columns are the H-Matrix scale factor which in this simulator 
1s 1.0. The remaining columns are left blank. The order of the runway 
points is as follows: far-centerline, far-left, near-left, near-center- 
line, near-right, far-right. 

Card 7 - The first ten columns consist of the optimum glide slope 
(GS0) in degrees. The second ten columns contain the x coordinate of 
the location of the Fresnel Lens (XMIR) in feet. The third ten columns 
contain the velocity of the headwind (WIND) in knots. 

Card 8 - 10 - These cards contain the inputed scale factors for the 
scaling of the analog solution. Each scale factor has a field of ten 
columns and are ordered as follows: 

P SP SFZ, STFS, SFA, 5-5, 5Ғ6, 5Ғ8, 559 

0 SEI] 5РІ2, shia, oFl4, SFIS, SF16, SF17 

Card 10 - SF18, SF19, SF30, the remaining columns are blank 

Card 11 - This card contains the initial conditions of the problem 
in ten column increments as follows: initial x coordinate (SX0) in feet, 
initial altitude (SZ0) in feet, initial velcoity (VO) in ft/sec, initial 
angle of attack (A10) in degrees, initial elevator angle (DITO) in 
degrees and the density of the air (RHO) in slug/ft°. The remaining 
portion of the card is left blank. 
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Card 12 - This card contains a portion of the aircraft constants 
again in ten column intervals. They are as follows: weight (W) of the 


aircraft in lb, wing span (B) in ft, wing chord (CB) in ft, wing area 
2 


3 


(S) in ft*, moment of inertia about the x axis (RIXX) in slug-ft 
moment of inertial about the y axis (RIYY) in slug-ft*, moment of inertia 
about the z axis (RIZZ) in slug-ft*, and the cross product of inertia 
(RIXZ) in slug-ft^. 

Card 13 - This card contains some more aircraft constants again in 
ten column intervals. They are as follows: maximum thrust (T) in lbs, 
optimum/landing angle of attack (AOALDG) in degrees, height of eye at 
touchdown (SZF) in feet. The remaining columns are left blank. 

Cards 14 - 85 (Alternating) - These cards contain the aerodynamic 
coefficients dependent on both alpha and beta. The first card in the 
series is as follows: in columns 1-2 the coefficient number, columns 
3-4 blank, column 5 the index on beta (IB), colums 6-10 blank. On the 
remaining portion of the card in ten column intervals are the coefficients 
indexed on alpha (IA). The second card in the series also contains the 
remaining coefficients again in ten column intervals and indexed on alpha. 
Alpha is indexed 15 times for each beta index while beta is indexed 9 
times for each of the four coefficients. The range on alpha is -25 
degrees to * 45 degrees with 5 degree intervals while the range on beta 
is -40 degrees to +40 degrees with 10 degree intervals. 


The coefficients are as follows: 


Number Coefficient 
] С, 
2 Cn 
3 Cn 
4 Cy 
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Cards 86 - 125 (Alternating) - These cards contain the aerodynamic 


coefficients dependent only on alpha. The first card in the series is 

as follows: columns 1-2 the coefficient number, columns 3-10 blank. The 
remaining portion of the card in ten column intervals are the aerodynamic 
coefficients indexed on alpha. The second card of the series contains 
the remaining coefficients again in ten column intervals and indexed on 
alpha. Alpha is indexed 15 times for each of the coefficients. A list 


of the remining coefficients and their numbers is given below. 


Number Coefficient 
5 Ср 
6 CL 
7 С 
E 
8 Cn 
Sit 
9 С 
Tor 
10 C 
lar 
11 C 
sr 
12 C 
m 
13 C 
Vea 
14 C 
ба 
15 С 
Үр 
16 С 
Тр 
17 "np 
18 C 
YR 
19 С 
In 
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Number Coefficient 


20 С 
IR 

2] С 
LQ 
22 C. 
Q 
23 C 
а 
24 C 
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APPENDIX G - SAMPLE OUTPUT 


The following is a sample output of the digital program and contains 


the following sections: 


Í. 


> jw N 


сл 


Runway Constants 

Scale Factors and Aircraft Constants 
Original Aircraft Aerodynamic Coefficients 
Usable Aircraft Aerodynamic Coefficients 
Pot Settings for the Analog Computer 


Output of the Amplifiers and D/A Trunks 
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T 
b 
pa 
T 


0% 


LCE SESE 
τε. 
Адый ү >л 


АРЕ ОА 


G Ge c 
OES *9 
CCO*OCEZ 


0.007 


c*o0t 


ΠΤΙ; 
NS Com 


(s 


А 


SINVISNOO 


= aM 
= EM 

з eM 

z Ти 

з «УМ 

= 659 

= IH 

= SIWX 

O° 9 
0.0008 S 
0.0008 + 
0*0008 E 
o € 
ο T 

X ЇХ Іба 
ДУМКИ 
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που ον 
020%09С001 
OOO. bcm 
ООО r5S/S9 
000 25285 
000%<4%2Т 
000%%5/6 
Orga 
DES 
000%58087 


οσο το” 
85000 * 

= 
076%" 
000*06d 
oCo OO I 
020526066 


оок 

81€ * 1 
Breet 
Ser 
986%0 
986%2 
090%! 
200%! 
Сео" 

Sc 

οσο 


- 


ZXIM 
Z271% 
ЛААЇЗ 
XS 


a3 


СМІМ 
QHA 
0110 
ОТУ 
Q^ 
O25 
OXS 


Od 
GE JS 
SESS 
Leas 
9685 
Sess 
8Е 35 
CES 
ΘΕ 39 
DES 
0€35 


009%0Т 
000%" 
OO η 
000%" 
009%Т 
900%Т 
οσο» ος 
ους. ος 
0002902 
000%1 
200° 90 
000 00" 
000*S 
000%< 
ОООО 
0060006 
00030007 
ο στο ος 
000%<8 
000%58 
0092%<8 
О0С%<" 
000558 
00076 
сое 
000 •2 
ООО е 
000100% 
000 00% 
000%00" 


3 


пи нини 


n ἡ 


αἼτεον 


5255 
Sols 
Leas 
9035 
50 35 
Ше 
Є с 35 
Ὁ ὁ 35 
tess 
0856 
6 7 35 
91315 
¿135 
91 39 
αἱ 35 
#145 
€ T3S 
2145 
1735 
Ο 135 
63S 
535 
{4% 
945 
С 35 
Nas 
E se 
Cas 
145 
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6644960" 
6645960" 
6655960" 
6655960" 
6655960" 
6654960" 
6655960" 
6655960" 
6655960" 
6695960" 
6655960" 
6655960" 
6655960"? 
6655960" 
6655960"? 


0° 0# 


66i#2L0° 
661*#2/{0* 
66T^24C* 
66l^2LC* 
66 Trato” 
66THZLO* 
661%240“ 
661240 * 
661 +240 * 
661%2¿0° 
66І%240% 
661%240% 
66іт220% 
66Tr2L0* 
66TH2L0* 


ODE 


00828%0% 
00982889“ 
00828,0” 
00828%0% 
0082889" 
00828480” 
00828+0* 
00828%0% 
00828802“ 
00828%0% 
00828%0% 
00828+0* 
00828+0* 
00828+0* 
00828%0% 


020g 


Cor T+20* 
004814920" 
00471720” 
0041920” 
00+TH20* 
оов тео“ 
00714720” 
οοπῖπσο» 
00474%20% 
00%7%20% 
00714720“ 
0041420% 
00%517%20% 
00%17%20” 
00+T9+20* 


0°07 


0200000* 
0000000* 
0000000* 
0000000• 
0000000 ° 
0000000* 
0000000* 
0000000% 
0000000“ 
0000000* 
0000000* 
0000000* 
0000000* 
0000000? 


0000000* 


0% 


T 830uNN 


904 TH20* 
00%1%420% 
0041420” 
0041720“ 
90%7720% 
00%Т7420% 
0041422: 
009 TH20* 
Q00+TAN20* 
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OO0*T*2C* 
O00*T*20* 
00%1%20% 
004 TA20* 


00%7%20% 


ο οι» 
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00828%0% 
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00828%0% 
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0082840: 
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00828%0» 
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661%240• 
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661+240* 
661%2¿0° 
66TH2LO* 
АЫ 
661+240• 
661%220% 
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661%с40• 
661+с40• 
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66TH2L0* 
66Th240° 


0%0Е- 
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6655960" 
66535960"? 
66559ь0" 
665596С" 
60655960" 
6655960"? 
66535960" 
66535960"? 
6655960" 
6655960" 
6655960"? 
6655960? 
6655960" 


6655960? 


0%0%- 


үналу 
VIS 


АЧЫЧУ 3392 IVNIEISG 
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